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energies. It is argued that the invariant mass dependence of HBT radii extracted from the correla- 
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quark gluon plasma expected to be formed in nuclear collisions at RHIC and LHC. Quantitatively 
different magnitudes of HBT radii are obtained at RHIC and LHC indicating stronger radial flow 
at LHC. 
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1. Introduction 

The nuclear collisions at ultra-relativistic energies have been effected at Relativistic Heavy Ion 
Collider (RHIC) and Large Hadron Collider (LHC) to create and study the properties of thermal 
phase of quarks and gluons which is called quark gluon plasma (QGP). Several observables have 
been proposed to study the properties of QGP [[T]] - among those the electromagnetic probes - 
photons and dileptons are advantageous because (i) they are produced at every space-time point of 
the QGP (ii) they leave the system without rescatterings and hence can deliver the information of 
the source point without any distortion ( for review). 

It is well know that the electromagnetic probes can be used as thermometer as well as a 
flowmeter too [Q, |5|]. The transverse momentum {kj) distribution of photons and dileptons reflect 
the temperature of the source. Because their productions from a thermal source depend on the 
temperature {T) of the bath through the thermal phase space factors of the participants of the re- 
action that produces the photons and dileptons. However, the thermal phase space factor may be 
changed by several factors - e.g. the transverse kick due to flow received by low kj photons from 
the low temperature hadronic phase will mingle with the high kj photons from the partonic phase, 
making the task of detecting and characterizing the QGP difficult. For dilepton the situation is, 
however, different because in this case we have two kinematic variables to describe the spectra - 
out of these two, the kj spectra of lepton pairs is affected but the kj integrated invariant mass (M) 
spectra is unaltered by the flow. Moreover, the M spectra of thermal dileptons is dominated by the 
late hadronic phase for M <mp and by the early QGP phase for M > m^, where nip and are 
p and <p mesons mass respectively. Therefore, the M spectra of dileptons can be used as a clock 
for heavy ion collisions - large M corresponds to early time and small M(~ mp) correlate to large 
time. This suggests that a judicious choice of kj and M windows will be very useful to map the 
temperature and flow of the evolving matter. 



2. Bose Einstein Correlation Function for Dileptons 

Experimental measurements of two-particle Hanbury-Brown Twiss (HBT) intensity interfer- 
ometry has been established as an useful tool to study the space-time evolution of the heavy-ion 
reaction [^]. The HBT interferometry with lepton pairs (or with virtual photons which differ form 
lepton pairs by a known factor) proceeds with the computation of the Bose-Einstein correlation 
function for two identical lepton pairs defined as, 

C2{kuk2) = -^S^ (2.1) 
m,)P,{k2) 

where kj is the three momentum of the pair, / and P\{ki) and P2{k\,k2) represent the one- and two- 
particle inclusive lepton pair spectra respectively and is expressed as follows: 



Pi{k) = / d'^x(oix,k) (2.2) 



and 



PiihM) = Pi{ki)Pi{k2) + - d'^xid'^X2 (o{xi,K)(o{x2,K) cos(Ajc''^^) (2.3) 
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where K = {k\+ ^2) /2, A/:^ = — ^2/1 = q^i, Ax^ = xi^ — X2^. x,-^ and are the four co- 
ordinates for position and momentum variables respectively. (o{x,k) is the source function related 
to the thermal emission of lepton pairs per unit four volume, expressed as follows: 

, , r^l 2 dR 
(0{x,k)= dM^—^-j—- (2.4) 

The inclusion of the spin of the lepton pairs (corresponding to the spin of virtual photon, which 
is 3) will reduce the value of C2 — 1 by 1/3. The correlation functions can be evaluated by using 
Eqs. U and U for different average mass windows, (M) ( —Mg+g-)= (Mi +M2)/2. 



The leading order process through which lepton pairs are produced in QGP is qq — )• /+/^ [^]. For 
the low M dilepton production from the hadronic phase the decays of the light vector mesons p , CO 
and ^ have been considered including the continuum [@,|8|]. Since the continuum part of the vector 
meson spectral functions are included, the processes like four pions annihilation [^] are excluded 
to avoid double counting. 



Table 1: Values of the various parameters used in the relativistic hydrodynamical calculations. 



Input 


RHIC 


LHC 


Nf 


2.5 


3 


Ti 


290 MeV 


640 MeV 




0.6 fm 


0.1 fm 


Tc 


175 MeV 


175 MeV 


Tch 


170 MeV 


170 MeV 


Tfo 


120 MeV 


120 MeV 


EoS 


2+1 Lattice QCD 


2+1 Lattice QCD 



For the space time evolution of the system relativistic hydrodynamical model with cylindrical 



symmetry [10] and boost invariance along the longitudinal direction [11] has been used. The 
initial temperature (7]) and proper thermalization time (t,) of the system is constrained by the 
hadronic multiplicity (dN/dy) through the relation dN/dy ~ Tf'rj. The equation of state (EoS) 
which governs the rate of expansion/cooling has been taken from the lattice QCD calculations [p^]. 
The chemical (Tch) and kinetic (Tfo) freeze-out temperatures are fixed by the particle ratios and the 
slope of the kj spectra of hadron [[l3|]. The values of these parameters are displayed in Table |I[ 

With all the ingredients mentioned above we evaluate the correlation function C2 for different 
invariant mass windows for 0-5% Au+Au collisions with 0-5% centrality at ^/s^n = 200 GeV and 
Pb+Pb collisions at ^/s^n - 2.76 TeV as a function of qside and q„ut which are related to transverse 
momentum of individual pair as follows (see [Q] for details): 

qout = {k\j -klj)/ f{kiT,k2T) 



q.ide = (IkirkiTX 1 -cos2(i/Ai - ^fi)) / f{kiT Mt) (2.5) 



where f{k[T,k2T) = \ + ^27 + 2^17^27 cos(i//'i — i//2)- 
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3. The HBT Radii 

The dimension of the source can be obtained by parameterizing the calculated correlation 
function (C2) of the lepton pairs by the empirical (Gaussian) form: 



C2 = i + -^M-RWi) 



(3.1) 



where the subscript / stand for out and side and A represents the degree of chaoticity of the source. 
The deviation of A from 1 will indicate the presence of non-chaotic sources. We have evaluated 
the C2 for (M) =0.3, 0.5, 0.7, 1.2, 1.6 and 2.5 GeV and the HBT radii are extracted by using 
Eq. 3J . The R^ide scales the transverse dimension and the Rout measures both the transverse size 



and duration of particle emission [14p. The extracted R^ide and R^ut (using Eq. |3.1| ) for different 



(M) are shown in Fig. [I] both for RHIC and LHC energies. 
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Figure 1: The variation of /?i,y(,(left panel) and 7?„„,(right panel) with (M) for RHIC (dashed line) and LHC 
(soUd line) energies. Rside is evaluated with k\j — kxr ~ 1 GeV and = and for y/i = y/2 = and 
kn = 1 GeV. 



In Fig |I] (left panel) we display the variation of Rside with M. It can be shown that the Rside is re- 
lated to the collective motion of the system through the relation: Rside ~ 1/(1 + ^collective /^Ihermal ) • 
With the onset of transverse expansion a rarefaction wave moves toward the center of the cylindri- 
cal geometry - as a consequence the transverse dimension of the emission zone reduces with time. 
Therefore, the size of the emission region is larger at early times and smaller at late times. The 
high (M) regions are dominated by the early partonic phase where the collective flow [||] is not 
developed fully i.e. the ratio of collective to thermal energy is small hence show larger R^ide of 
the source. The lepton pairs with M mp are emitted from the late hadronic phase where the size 
of the emission zone is smaller due to larger collective flow effects giving rise to a smaller R^ide- 
The ratio, Ecoiiective / Ethemai is quite large for M ^ nip and hence Rside is small which is reflected 
as a dip in the variation of R^ide with (M) around the p-mass region (Fig. |l| left panel). Thus the 
variation of R^ide with M can be used as a yardstick to scale the size of the evolving matter. The 
change of Rside with (M) for RHIC and LHC is qualitatively similar but quantitatively different. 
The smaller values of Rside for LHC is due to the larger radial expansion which can be understood 
from the fact that the quantity Ecoiiective / Ethermai is larger at LHC than RHIC. 
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The Rout probes both the transverse size and the duration of emission. The large M regions 
are populated by lepton pairs from the partonic phase where the effect of flow is low (size is 
large) but the duration of emission is small - resulting in a small values of Rout- Lepton pair with 
M ^ nip suffers from larger flow effects which should have resulted in a minimum as in /?s,yg in 
this M region. However, Rout probes the duration of emission too, which is large for hadronic 
phase because of the slower expansion due to softer EoS used in the present work for the hadronic 
phase. The larger duration compensates the reduction of Rout due to flow resulting is a bump in 
Rout for M (Fig. |T]right panel). The Rout at LHC is smaller than RHIC because the larger flow 
(corresponds to smaller size) at LHC compensates other factor (like duration of emission) which 
has an enhancing effect on Rout- 

The Rout and /?s,ye are the measures of the average size of the system [?] and has some de- 
pendence on space-time evolution models that is used. However, in the ratio. Rout /Rude some of 
the uncertainties associated with this modeling get canceled out. The ratio, Rout/R&ide gives the 



duration of particle emission [ ]15| , |16| , |17[ ] for various domains of (M) corresponding to different 
time slices. 




Figure 2: The ratio RomlRskie (left panel) and the difference \J R^^t — R\de ^^^^^ panel) as a function of (M) 
for RHIC (dashed line) and LHC (solid line) energies is shown. 



Fig. ^ depicts change of the RoutlRside and Rdiff{= \jR^ut ~R^nde) ^ function of (M) for 
RHIC as well as LHC energies. Both show a non-monotonic dependence on (M). The results 
reveal that both the ratio and the difference of HBT radii for LHC at low M are larger than the 
corresponding quantities at RHIC. 

Now we indicate below the experimental challenges in performing the HBT interferometry 
with lepton pairs at LHC (for RHIC see We compute the number of events from the luminosity 
(Jf ), the pp in-elastic cross section (a) and the run time {^) of the LHC as: 



50 X 10^^/(cm^.sec) and a=60 mb we get A'^ei 



(3.2) 

2 X 10^°. As an 



For j7 = 12 weeks, ^ 

example for (M) =500 MeV and kj = 1 GeV, the value of (dN/d^hrdy) for Pb-i-Pb coUision at 
^/5Kw=2.76 TeV is ~ 0. 138 x 10"^ Therefore, the total (for 2 x 10^° number of events) differential 
number of pairs in the above range of kj and M is ~ 2 x 10'" X 0.138 xlO-^~ 2.7 xlO^ 
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Similarly for the (M) =1.02 GeV and kT=l GeV The total (differential) number of pairs is 
2 X 10^. In this domain of kj and M the number of pairs produced per event is ~ 10~^. Therefore, 
the probability to get two pairs is 10^^, Therefore, roughly 10^ events will be necessary to perform 
the interferometry with lepton pairs in this region of Ut and M. 

4. Summary 

In summary the correlation functions for dilepton pairs has been evaluated and the HBT radii 
have been extracted for both Au+Au collisions at RHIC and Pb+Pb collision at LHC energies for 
different (M) windows. We argue that the variation of HBT radii with M for dilepton pairs can be 
used as an efficient tool to follow the change of the dimension of the evolving system with time. 
The quantitative difference in the HBT radii at RHIC and LHC indicate the development of larger 
flow at LHC compared to RHIC. 
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